All-order evaluation of weak measurements 

— The cases of an operator A which satisfies the property — 1 — 
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Some exact formulae of the expectation values and probability densities in a weak measurement 
for an operator A which satisfies the property A^ = 1 are derived. These formulae include all- 
order effects of the unitary evolution due to the von-Neumann interaction. These are valid not 
only in the weak measurement regime but also in the strong measurement regime and tell us the 
connection between these two regime. Using these formulae, arguments of the optimization of the 
signal amplification and the signal to noise ratio are developed in two typical experimental setups. 
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I. INTRODUCTION 

Since the proposal of the weak measurement by 
Aharonov, Albert, and Vaidman (AAV) [l| in 1988, 
weak measurements have been investigated by many re- 
searchers. The idea of weak measurement has been used 
to resolve fundamental paradoxes in quantum mechanics 
such as Hardy's paradox [2] . In addition to many theoret- 
ical works on weak measurements, it is important to note 
that some experiments realized this weak measurement in 
different experimental setups Q . These experiments 
show that the weak measurement is also very useful for 
high-precision measurements. For example, Hosten and 
Kwiat used the weak measurement to observe a tiny 
spin Hall effect in light; Dixon et al. d 13 (DS JH) used 
the weak measurement to detect very small transverse 
beam deflections. The original AAV work [Ij] also in- 
cludes the proposal of the application to the sequence 
of the Stern-Gerlach experiments for spin-1/2 particles. 
They claim that we can observe the spin of particles as 
a larger value than the range of its eigenvalues. This 
is called "weak- value amplification". The above high- 
precision measurements using the weak measurement are 
due to the effect of this weak- value amplification. 

Weak measurements are based on von-Neumann's 
measurement theory Q in which the total system con- 
sists of the system to be measured and a detector to 
measure the system. Further, we specify the initial state 
(pre-selection) and the final state (post-selection) of the 
system. AAV also proposed the situation of the measure- 
ment, in which the initial variance in the momentum con- 
jugate to the pointer variable of the detector is so small 
that the interaction between the system and the detector 
is very weak [To]. Because of this weakness, the measure- 
ment proposed by AAV is called "weak measurement". 
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In the linear- order of the interaction between the sys- 
tem and detector, the outcome of the weak measurements 
is so-called "weak value" . The weak- value amplification 
is essentially due to the fact that the weak value of an 
quantum observable may become larger than eigenvalues 
of this observable when the pre- and the post-selection is 
nearly orthogonal. Due to the weakness of the interaction 
between the system and the detector, the measurement 
by a single ensemble is iniprecise. However, as noted by 
Aharonov and Vaidman [10], the measurement become 
precise by a factor ^/N through performing large N en- 
semble experiments. 

Measurements of arbitrary strength beyond the linear- 
order interaction has been first discussed by Aharonov 
and Botero [ll| in the context of the framework called 
"Quantum average of weak value." In this framework, 
the strong measurement of a pointer variable can be 
regarded as quantum superpositions of weak measure- 
ments. They applied their framework to a specific case 
of a spin measurement. Furthermore, all-order effects of 
the unitary evolution due to the von-Neumann interac- 
tion between the system and the detector are also in- 
vestigated by investigated by Di Lorenzo and Euges [l^ 
in AAV setup to clarify the detector dynamics in weak 
masurements. 

More recently, Wu and Li [l^ proposed the general 
formulation of the weak measurement which includes all- 
order effects of the unitary evolution due to the von- 
Neumann interaction between the system and the detec- 
tor. Through this formulation, they took some higher- 
order effects into account when they computed the shift 
of pointer variables and pointed out that there is a over- 
lap of the pre- and the post-selection at which the out- 
come of the weak measurement have the maximal am- 
plification. However, since they did not take all higher- 
order effects into account, their claim on the maximal 
amplification is weak. 

In this paper, we carry out the all-order evaluation 
of some expectation values of pointer variables after the 
post-selection based on the formulation proposed by Wu 
and Li [l^ . Although the all-order evaluations of the ex- 
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pectation values in general weak measurement are diffi- 
cult, these evaluations are possible if we concentrate only 
on the weak measurements for an operator A of the sys- 
tem which satisfies the property = 1. Choosing the 
initial state of the detector as a zero mean- value Gaussian 
state, we derive some formulae of the expectation values 
and probability densities for the detector after the post- 
selection without any approximation. Through these for- 
mulae, we discuss the maximal amplification which sug- 
gested by Wu and Li. 

Although our consideration is restricted only to the 
case of the weak measurement for an operator A which 
satisfies the property A^ = 1, this case includes many 
experimental setups. For example, the weak measure- 
ment of the spins of spin- 1/2 particles, which was orig- 
inally proposed by AAV is included since the Pauli 
spin matrices satisfy the property A^ = 1. The experi- 
ment by Hosten-Kwiat Q and the experimental setup by 
DSJH are also included in our case, though there are 
some additional modification in their actual experimen- 
tal setups. Thus, our consideration will be applicable to 
many experimental setups. Therefore, it is worthwhile to 
research the weak measurements for an operator A which 
satisfies the property A^ = 1. 

Furthermore, we note that some experiments of weak 
measurement for an operator A which satisfies the prop- 
erty A^ = 1 are classified into two types: one is the weak 
measurements with a real weak value; and the other is 
those with a weak value of pure imaginary. A typical ex- 
ample of the weak measurement with a real weak value 
is the experimental setup of a spin- 1/2 particle proposed 
by AAV [IJ. On the other hand, a typical example of 
the weak measurement with a weak value of pure imag- 
inary is the DSJH experiment @. We apply our results 
of all-order evaluations to these two specific experimen- 
tal setups. Then, we discuss the optimizations of the 
expectation value of the pointer variable of the detector 
(i.e., the signal optimization) and the optimization of the 
signal to noise ratio (SNR). Through these applications, 
we concretely discuss the maximum amplification in the 
weak measurements. 

Organization of this paper is as follows: In Sec. |lTl we 
briefly review the general formulation proposed by Wu 
and Li. In Sec. IIIIl we summarize the formulae for some 
expectation values and probability densities which are 
derived from all-order evaluations through Wu-Li formu- 
lation. In Sec. lIVi the application of our formulae to AAV 
setup is discussed. In Sec. |Vl we discuss the application 
of our formulae to DSJH setup, though the experimental 
setup in this paper is a simpler version of the original 
DSJH setup. Final section (Sec. IVip is devoted to the 
summary. 

Throughout this paper, we use the natural unit h = 1. 



II. WU-LI FORMALISM 

Here, we review the description of weak measurements 
proposed by Wu and Li [isj. In Sec. Ill Al we first review 
the general framework of the weak measurement follow- 
ing Ref. [ll]. To carry out the analyses, we must treat 
two cases separately for a technical reason. One is the 
case where the initial and the final states of the system 
is not orthogonal, which is described in Sec. IIIBl The 
other is the case where the initial and the final states of 
the system is orthogonal, which is described in Sec lIICI 

A. General framework 

The total system we consider here is described by the 
density matrix p = ps ® Pd- is the density matrix of 
the "system" which is a quantum system and we measure 
an observable A associated with this system, pd is the 
density matrix of the "detector" which interacts with the 
system through the von-Neumann interaction 

gS{t^to)A(g>p. (2.1) 

Here, p is the conjugate momentum to the pointer 
variable q of the detector, i.e., [q,p] = i. In the 
usual von-Neumann interaction (strong interaction), the 
eigenvalues of A appear in the pointer variable q 
Using these three elements, the weak measurement 
is carried out through the sequence of four measure- 
ments p. First three processes of these four mea- 
surements are called "pre-selection" , "weak interaction" , 
"post-selection" . The final one is the measurement of the 
detector pointer variable through any type of the mea- 
surement in quantum mechanics. 

First, we prepare the initial state ps of the system 
through the projection measurement at t < to, which is 
called "pre-selection". We also prepare the initial state 
of the detector pd- After this pre-selection, the system 
and the detector interact with each other through the in- 
teraction Hamiltonian (12. ip . The time evolution through 
this interaction is described by the evolution operator 
U = e^^^^'P and the total density matrix p evolves as 

p' ^ UpU^ =P+Y ^^^ad>A o p, (2.2) 

n=l 

where ad" for arbitrary operators and O is recursively 
defined as 

ad^noQ := adf? o 9 = [O, 6] , (2.3) 
adrrioQ := adf? o (ad"-^17 o 9) 

= [f7,ad""il7oe] . (2.4) 

The prime in Eq. (j2.2p denotes the operator after the 
interaction (|2.ip . The density matrix of the system after 
this interaction is given by 

= Tr^p' = Ps + J2 ^^^(p")ad"A o p„ (2.5) 
^ — ^ nl 

n—l 
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where Tr^ means taking the trace of the detector den- 
sity matrix and (p") :— Tid{p^Pd)- Equation (|2.5p 
imphes that the density matrix of the system hardly 
changes through the interaction with the detector if 
(7sup{(p")^/",n G N} ^ 1. Roughly speaking, this 
condition is regarded as gAp <C 1, where (Ap)^ is the 
variance in p, and is interpreted that the interaction be- 
tween the system and the detector in the measurement 
is "weak interaction" . After this interaction, we restrict 
the final state of the system by the projection opera- 
— > /o'llf. This restriction is called "post- 



by 



tor n 



selection" . The density matrix of the detector after the 
post-selection is given by 



Pd — 



Trp'H/ ' 



(2.6) 



where Tr^ (Tr) means taking the trace of the system den- 
sity matrix (the total density matrix). 

Although ps and pd may describe mixed states of the 
system and the detector, we restrict our attention to pure 
states as the initial density matrices ps and pd- We de- 
note these initial density matrices as ps = ^-nd 
Pd = If/*) (</'!■ Further, we also denote the projection op- 
erator for the post-selection by 11/ := |-0/)(?/'/|. In this 
case, the normalization factor (Trp'II/) of the density 
matrix of the detector after the post-selection is given by 



+ 00 



n=l k=0 

x(V^/|A"-'=|V'.)(V'.|A'=|V/),(2.7) 



where (p") = ((/)|p"|(^) and „Cfc is the binomial coeffi- 
cient. 

To carry out the further analyses, the factor 
plays an important role and separate treatments are re- 
quired according to the fact whether |(V'/|V'i)| = or 
not. 



Non-orthogonal weak measurement \{tpf\i}}i)\ 7^ 



Here, we consider the case where the pre- and post- 
selection are not orthogonal, i.e., \{ipf\^i}\ 7^ 0. In this 
case, the trace of the post-selected density matrix and 
the density matrix p'^ after the post-selection are given 



+ 00 



(2.8) 



Harip") 



^(-l)^,Cfe(A"-'=)„(A'=) 



fe=0 



(2.9) 



+ 00 



n=l k=0 

x(A"-^-).^(A''^)>"-Vdp', 

(2.10) 

where (•)„ := {ipf\ ■ 

When the wave function (p|(/)) is even inp, i.e., (p") — 
for odd n, Wu and Li derived the formulae of the shifts 
in q and p as 

,g5RA„ +g3A^({g,p}) 
= [ — — z : — T' (2-11) 



l + g^p^){\A^f-n{A^) 



2,g3A,„(p2\ 



Sp 



l + g^P^) |A^|2-5R(A2 



(2.12) 



where Sq TT{qp'^) — Tr{qpd) and Sp :— Tt{pp'^) — 
TT{ppd). In their derivation, they neglect terms of 0{g^) 
in the numerators and the denominators, but they do not 
expand the total expressions (|2.1ip and (|2.12p in form of 
the power series of g. Although these treatments of Sq 
and Sp might be regarded as some renormalization tech- 
nique, it is also true that the expressions (|2.1ip and ()2.12p 
include only partial effects of higher order of g. 

As pointed out by AAV [l|, weak values may become 
very large in the limit \{fpf\tpi) \ — >■ (7^ 0). At the order 
of 0{g), the shifts (12.111) and (I2.12p are proportional to 
the weak value (A)^, [IJ]- This implies that the shifts 
(12. lip and (12.121) of order 0{g) may diverge in the limit 
KV'/lV'i)! ~^ 0. This is the essence of the weak value 
amplification. However, from the total expressions of 
Eqs. (|2.1ip and (|2.12p . Wu and Li suggested that, in the 
limit |(V'/|'0i)l ~^ 0, these shifts decrease rapidly when 
|Au,|^ become comparable with [g"^ ip^))^^ ■ This argu- 
ments implies that, for a fixed 5^(p^), there may exist 
a maximum shift of a pointer quantity, and an optimal 
overlap |('0/l'0i)l to achieve the maximum shift. We call 
this overlap as the optimal pre-selection (or optimal post- 
selection). 

Although Wu and Li claim is weak in the sense that 
they did not take all higher-order effects into account, in 
this paper, we show that their claim on the optimal pre- 
selection is essentially correct through the all-order eval- 
uation of weak measurements for an operator A which 
satisfies the property A^ = 1. 
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C. Orthogonal weak measurement ~ 



A. Non-orthogonal weak measurement litpflipi}] 7^ 



Next, we consider the orthogonal case where 
= 0- In this case, the original formalism of 
the weak measurement fails and the weak values are not 
defined. This is easily seen from the fact that the nor- 
malization factor Z defined by Eq. (|2.8p is ill-defined. 
However, instead of Z, Wu and Li defined Zo by 

TTp'Uf=:g'{p')\{^f\A\^,)fZ,, (2.13) 



° " n! 

n— 1 ^ ' 



fc=0 



where 



(^/|A"+i|^, 



(2.14) 
(2.15) 



Wu and Li called (A") on, defined by Eq. (|2.15l) as or- 
thogonal weak values. The density matrix of the detector 
after the post-selection is given by 

Zo{p^)p'd = PPdP 

n=l ■ k=0 

xp"-"+V<i/+'- (2.16) 

From this expression (I2.16p . Wu and Li claim that the 
orthogonal weak values (I2.15P play the similar role to 
the original weak values in non-orthogonal case. 



III. ALL-ORDER EVALUATION OF WEAK 
MEASUREMENTS FOR AN OPERATOR A 
WHICH SATISFIES A^ = 1 

Here, we evaluate the density matrix of the detector 
after the post-selection and some expectation values in 
the case for an operator A which satisfies the property 
A^ = 1 based on the Wu-Li formalism. In addition to 
the restriction of our consideration to the simple operator 
case, in this section, we assume that the initial state pd = 
|(/))(^| of the detector is zero mean-value Gaussian, i.e.. 



1 



27r(p2 



1/4 



exp 



4(p2 



(3.1) 



From this initial state of the detector, we can easily derive 
the properties of the initial state: 

(p2"+i)=0, (p2") = (2n-l)!!(p2)". (3.2) 

As reviewed in the last section [TTl according to the 
norm have to treat the density matrix in 

different way. Therefore, we treat a non-orthogonal weak 
measurement and an orthogonal one, separately. 



When the initial state of the detector is zero mean- 
value Gaussian p.ip . the moments of p are given by 
Eqs. p.2p . In this case, the normalization Z [Eq. (jAll) ] 
is given by 



Z = 1 + 1(1-|(A) 



where s is a parameter defined by 



-1 



(3.3) 



(3.4) 



Similar calculations lead the expectation values of p and 
q after the post-selection 



9 

9{P)' 



Z ' 
se-"3(A) 



(3.5) 
(3.6) 



[Here, we denotes the expectation value of * for the 
detector after the post-selection by (*)'. Fluctuations 

Ag := ^{{q^{q)'f)' and Ap := ^W^WW in P 
and q after the post-selection are given by 



(Ag)2 



i + 2^(i + i<^>-i^ 



(3?(A)^)' 



Z2 



2Z 



l-l(A) 



(3(A).)2 



(3.7) 



(3.8) 



Further, the probability densities in p-space and q- 
space are given by 



{P\P» 

Wdk) 



l-l(A)^ 
25(A)^sin(2gp) 

2 



cos{2gp) - 1) 

{p\Pd\p} 
2Z ' 



(3.9) 



1-|(A 

1 + 



(A)„,Pl cosh ' ^'^'^ 



-H25R(A)„sinh (jy^ ^ 

^e~-H/l\pd\q) 
2Z 



(3.10) 



where {p\pd\p) and {q\pd\q) are Gaussian initial probabil- 
ity densities 



(p\Pd\p) 
{<l\Pd\q) 



■ exp 



(gp)' 



■ exp 



(3.11) 
(3.12) 
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The derivation of Eq. (I3.10|) is explained in Appendix 

ixn 

Here, we note that the parameter s defined in Eq. p.4|) 
is a measure of the strength of the interaction. Usually, 
it is said that the interaction between the system and the 
detector is weak if the coupling constant g is very small. 
On the other hand, in the weak measurement ||10], it 
is said that the interaction between the system and the 
detector is weak if the initial variance of the pointer vari- 
able q is very large, i.e., the variance in the conjugate 
momentum p is very small. These two concepts of the 
"weakness" of the measurement are automatically repre- 
sented by the single non-dimensional parameter s. We 
call s as the coupling parameter, and say that the inter- 
action between the system and the detector is weak if 
s <C 1 and strong if s 1. 

We have to emphasize that our formulae shown here 
are the results from the all-order evaluation of s and valid 
not only in the weak measurement regime s ^ 1 but also 
in the strong measurement regime s 3> 1. The results 
coincide with those of the measurement in the strong 
regime. This situation can be observed through the spe- 
cific experimental setups discussed in Sec. IIVI 

Finally, we note that probability distributions for weak 
measurements (both in the strong and weak regime) are 
first discussed by Aharonov and Botero flTl in the con- 
text of the framework called "Quantum averages of weak 
value" as mentioned in Sec. H] Of course, our formulae of 
the probability distribution shown in this paper are not 
general because we concentrate only on the case of an 
operator A which satisfies the property A^ — 1. How- 
ever, as emphasize in Sec. HI many experimental setups 
are included in this special case and we have derived ex- 
plicit simple analytic formulae for this special case. This 
is one of main points of this paper. 



B. Orthogonal weak measurement |(V'/IV'j)l — 

Now, we consider the orthogonal case where the pre- 
selected state and the post-selected state are orthogonal 
to each other, i.e., {ipf\ipi) = 0. As reviewed in Sec. IH Cl 
the density matrix of the detector after the post-selection 
is given by Eq. (j2.16p . In the case of the weak measure- 
ments for the operator A with the property A^ = 1, the 
orthogonal weak values (|2.15p are given by 



1 



for n is even, 
n + 1 
for n is odd. 



(3.13) 



This expression implies that no information of A appears 
in the orthogonal weak measurement for an operator A 
with the property A^ = 1. 

Through the Gaussian initial state p.ip of the detector 
with the properties p.2p . the normalization constant Zg 
defined by Eq. ()2.14p is given by 



(3.14) 



where the coupling parameter s is defined in Eq. (13.41) . 
The density matrix (|2.16p of the detector after the post- 
selection is given by 



- (2n- 



fc=0 



(3.15) 



In the case where the initial state of the detector is 
zero mean- value Gaussian p.ip . the expectation value of 
p after the post-selection, which is evaluated in Appendix 
21 is trivial. 



{pY = 0, 



(3.16) 



due to the properties p.2p . Further, in Appendix lA 21 
we also show that the expectation value of q [Eq. (jA33p ] 
after the post-selection also yields a trivial result 



(g)' = 0. 



(3.17) 



As shown in Appendix lA 21 the fluctuations Ap and 
Aq in p and q of the detector after the post-selection are 
given by 



9\^pf = 
{^qf 



s 1 



(25 - l)e- 



2 4 - 3s - Ae-'' 
1 1 - g-'' + 4s 
2s4-4e-^ -3s' 



(3.18) 
(3.19) 



The probability densities in p-space and in g-space are 
given by 



{P\p'a\p) 

(9|Pdk> 



1 — cos(2(7p) 
2 (4 - 46-^* - 3s) 
2e~^ sinh^ (sq/g) 



4-4e- 



3s 



{p\pd\p), (3.20) 
{q\pd\q), (3.21) 



respectively. Here, {p\pd\p) and (q\pd\q) are the Gaus- 
sian initial probability densities p. lip and p.l2p . re- 
spectively. The derivations of these formulae are given 
in Appendix I A 21 

Thus, both in the non-orthogonal weak measurements 
(Sec. nil Ap and the orthogonal one (Sec. IIII Bp . we ex- 
plicitly derived the analytical expressions of the expecta- 
tion values of p and q, fluctuations in p and q, and the 
probability distributions in p-space and in g-space for the 
detector only under two assumptions, i.e., the operator 
A for the system satisfies the property A^ = 1 and the 
initial state of the detector is zero mean-value Gaussian 

We note that the formulae p.Sp and p.6p for the ex- 
pectation values for p and q coincide with Eqs. (|2.1ip and 
(|2.12p . respectively, if we ignore the higher-order terms 
of than 0{g^). In this sense, equations (|3.5p and (|3.6p 
are all-order extension of Eqs. (j2.1ip and (I2.12p derived 
by Wu and Li fl^. We also note that the expressions 
of Eqs. p.Sp and (13. 6p are valid for arbitrary value of 



6 



the coupling parameter s. Furthermore, the behaviors 
of Eqs. p.Sp and p.6p are quahtatively same as those of 
Eqs. ()2.11|) and (|2.12|) . Therefore, we may say that the 
claim on the optimal post-selection proposed by Wu and 
Li is essentially correct and mathematically justified by 
Eqs. ([331) and ([3^]) . 

In the following two sections, we apply the formulae 
summarized in this section to two specific experimental 
setups and examine the weak measurement of these two 
setups in detail. Since we already showed that the or- 
thogonal weak measurement yields trivial results in the 
expectation value oi p and q, we concentrate only on the 
non-orthogonal weak measurement. 



IV. APPLICATION TO AAV SETUP 

In this section, we apply our formulae derived in 
Sec, mil to the AAV ^1,] setup. Through this application, 
we discuss the optimization of the expectation value of 
the signal and SNR. 



A. Setup of experiment 

The experimental setup proposed by AAV [Tj is the 
sequence of three Stern-Gerlach experiments for spin-1/2 
particles. 

The pre-selected state of the spin-1/2 particle is 

I tc> = ^ (Vl + sina| t.) + Vl-sina| , (4.1) 

which is an eigenstate t^) ~ +| t^) of the operator 
CTj = cos aax + sin aaz ■ 

The weak interaction in the weak measurement is de- 
scribed by the interaction Hamiltonian 



n 



-gz<g)azS{t-to), g ^ IJ. 



dz 



(4.2) 



where /i is the magnetic moment of the spin-1/2 particle 
and Bz is the z-component of the magnetic field. The 
pointer variable in this setup is pz which conjugate to 
z. We note that the operator A to be observed is the 
spin z-component of a spin-1/2 particle through the von- 
Neumann interaction (j2.ip . i.e., 



A cr^ 



t.)(t. i-u.)(;. 



(4.3) 



which satisfies the property A^ = 1. Then, we may apply 
our formulae provided in Sec. lIII Al Comparing Eq. (14. 2p 
with the interaction Hamiltonian (j2.1l) . we find the cor- 
respondence of variables as 



p-^~z, q^Pz- 
The post-selection in this setup is 

lt.) = ^(lt.) + U.)) 



which is an eigenstate (Jx\ tx) = +| tx) of the x- 
component ax of the spin. 

The weak value in this setup is given by 



(t.lt,) -*""2' 



(4.6) 



where a is the pre-selection angle. We note that this 
weak value (14.61) is real. 



B. All-order expectation values and probability 
distribution 



Here, we apply the formulae summarized in Sec. IIII Al 
to the AAV setup. The normalization factor Z 
[Eq. p.3|) ] , is given by 



cos a 



1 



cos a 



(4.7) 



where s is the coupling parameter [see Eq. (|3.4|) ] defined 

by 



Expectation value of pz and z are given by 



sm a 



1 



' cos a 



{zY = 0. 



(4.8) 



(4.9) 



The expectation value of Ipz)' / g is shown as a function of 
the coupling parameter s and the pre-selection angle a in 
Fig. [TJ In Fig. [DJa), we can see a pole at (s, a) — (0, tt). 
Due to this pole, the weak value is amplified as pointed 
out by AAV. In the region s ^ 1, the expectation value 
(|4.9p of {pzY /g behaves {pzY /g ~ sin a. This behavior 
can be seen in Fig. [TJc). The qualitative difference be- 
tween Eqs. dHH-dHH) by Wu-Li and Eqs. in 
this paper becomes large in the strong region s 3> 1. 
Fluctuations Ap^ and Az are given by 



(Ap.)' 
9^ 



1 cos a (cos a -|- e ") 



2s 



(1 -t- e *cosa)^ 



s s^e *cosa 
2 1 + e^" cos a 



, (4.10) 



(4.11) 



We also note that the first term in Eq. (liTTO)) [Eq. (|4?TTt ] 
shows the initial variance in pz (in z). The remain- 
ing terms in Eqs. (|4.10p and (|4.1ip are due to the pre- 
selection, weak interaction, and the post-selection. 

The probability density of the detector after the post- 
selection in p^-space is given by 



(4-4) {Pz\p',\Pz 



(4.5) 



cosa -I- cos\i{2spz/ g) + sinasinh(2sp2/5) 
1 -|- e"** cos a 



exp[-s] J;^exp [-s{pz/gf] 



(4.12) 
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FIG. 1: {Pz)' /g is shown as a function of the coupling s — 2g'^{z'^) and the pre-selection angle a in various range of s. (a) 
s £ (0,0.1), (b) s £ (0, 1), (c) s £ (0, 10). From these figures, we can see that, for given coupling parameter s, there is a optimal 
pre-selection angle a such that the expectation value of pz is maximized. 



Lorenzo and Egues [1^ also derived analytical formu- 
lae of the expectation of the pointer variable and the 
probability distribution in more complicated form. Their 
derivation is based on Horn's rule of the joint probability. 
Our results shown here are consistent with their results. 



C. Expectation value optimization 

Figure [T] explicitly shows the existence of the ridge in 
the surface of the expectation value (j4.9p . This means 
that for a fixed coupling parameter s, there is the optimal 
pre-selection angle a at which the expectation value (|4.9I) 
is maximized. This was pointed out by Wu and Li jl3l | 
from the less accurate expression (|2.1ip . On the other 
hand, we can accurately discuss this optimization of the 
expectation value from our exact expression (|4.9p . Here, 
we consider this optimization of the expectation value 
(1121) in detail. 

To derive the points at which the expectation value is 
optimized, we consider the equation d{{pz)' / g) / da = 0. 
This equation yields 



(4.13) 



We call the line which is expressed by Eq. (|4.13p on the 
(s,a)-plane as the optimal expectation- value line. On 
this optimal line, the expectation value of Pz and the 
fluctuation Ap^ in pz are given by 



9 

^Pz 

9 



1 



1 



^27' 



(4.14) 
(4.15) 



Here, we note that the fluctuation Ap^ (|4.15p in pz co- 
incides with that for the initial state of the detector. 
The optimal expectation- value line, the expectation value 
(pz)', and the fluctuation Apz in pz on this optimal line 
are shown in Fig. [51 



The expectation value (|4.14p of Pz explicitly shows that 
we can accomplish the arbitrary large weak value am- 
plification if we prepare the sufficiently small coupling 
parameter s. Actually, when s <C 1, {PzY ~ 9l^/^ — 
l/{2^/{z^). Thus, the expectation value of can be 
very large if we choose the initial variance in z is very 
small. This is just the weak value amplification proposed 
by AAV i. 

On the optimal expectation- value line (|4.13p . the prob- 
ability distribution (|4.12p is given by 



{Pz\p'd\Pz) = 



cosh(2sa;) + Vl — e "^^ sinh(2sa;) — e 



1 -e 



-2s 



(4.16) 



where x := Pz/g- This probability distribution on the 
optimal expectation-value line is shown in Fig. [3l with 
some coupling parameters s. In Fig. |3l s = 0.1 case 
corresponds to the weak measurement on the optimal 
expectation- value line (|4.13p . This shows that the peak 
of the probability distribution slightly deviates from the 
weak value A^^ and the probability density after the 
post-selection is slightly different from the Gaussian pro- 
file [HI. s = 1.0 case is still essentially same as s — 0.1 
case, s = 10 and s = 1000 cases correspond to the strong 
measurement regime. On the optimal expectation-value 
line (|4.13l) . a approaches to 7r/2 in the limit s — >■ cx) 
and Apz in Eq. (|4.15l) approaches to 0. Here, we note 
that the pre-selected state with a = 7r/2 corresponds 
to the eigenstate | fz) of Uz- In this case, we measure 
this eigenvalue -|-1 with small uncertainty. This situa- 
tion is well-described by the behavior of the probability 
distribution with s — 1000 in Fig. [31 Therefore, the prob- 
ability distribution (|4.16p well-describes not only in the 
weak measurement regime s ^ 1 but also in the strong 
measurement regime s 3> 1. 

Although we have an arbitrary large expectation value 
(|4.14p if we choose s <C 1, small coupling parameter 
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FIG. 2: [Online Color] The optimal expectation-value line 
(|4.13p [top panel], the expectation value (pz)' /g [Eq. (|4.14[) ] 
and the fluctuation Ap^ [Eq. (|4.15l) ] in pz after the post- 
selection on the optimal line (|4.13ll [lower panel] are shown 
as functions of the coupling parameter s. We also show 
that ((pz)' ± l^Pz)/g in this lower panel. [The red line 
is ((pz)' + /^Pz)/g and the blue hue is ((pz)' — Apz)/(/.] 
This figure shows that if we choose the small parameter of 
s ^ 1, we can accomplish the large expectation value of 
(pz)'/5 ~ l/\/2s, but the fluctuations in pz also amplified 
as Apz/g ~ 1/a/2s. 



s gives large variance in pz, as shown in Eq. (j4.15l) . 
Actually, fluctuation Ap^ in Eq. (14.15^ also behaves as 
Apz = gl\f2s = l/(2-\/(i^). Since the fluctuation Ap^ 
is regarded as a noise in this weak measurement, this 
means that the SNR on the optimal expectation-value 
line (|4.13p is {p^)' / ^Pz ~ 1- Thus, we do not have a 
large SNR in the expectation-value (signal) optimization 
of the single particle experiment. Therefore, we consider 
the optimization of the SNR in the next subsection. 



D. SNR optimization 

The expectation value {pz)' [Eq. (14. 9p ] and the fluc- 
tuation Ap^ := ^{{pz - (pzYY)' [Eq. dUni)] after the 
post-selection are regarded as the signal and a noise in 
the measurement of pz- Therefore, in this section, we 



9^z\p'd\p^) 





FIG. 3: [Online Color] Probability distribution functions 
(|4.16p in Pz/ g-spa.ce on the optimal expectation- value line 
(|4.13p with some coupling parameters s are shown by the 
thick lines (red lines). We also plot the initial probability 

distribution y^e"'''''^''^-' with the same coupfing parameter 
s by the thin lines (blue lines), s — 0.1 case corresponds to 
the weak measurement on the optimal expectation-value line 
(|4.13p . This shows that the peak of the probability distri- 
bution slightly deviates from the weak value A^- s — 1.0 
case is still essentially same as s = 0.1 case, s = 10 and 
s — 1000 cases correspond to the strong measurement case. 
These behaviors of the probability distribution (|4.16|l also 
well-describes the strong measurement regime s S> 1. 




FIG. 4: The signal to noise ratio (SNR) (|4.17p is shown as 
a function of the coupling parameter s and the pre-selection 
angle a : (a) s G (0,1) and a € {n/2,n); (b) s G (0,10) 
and a £ (7r,0). There is the ridge of the SNR from (s, a) = 
(0, tt) to (s. Of) = (-1-00, 7r/2). [We note that the direction 
of s-axes are opposite to those in Fig. [Tj] The ridge around 
a = 7r/2 is due to the strong measurement regime s 3> 1, 
which is due to the fact that a = n/2 corresponds to the 
eigenstate | fz) of the operator CTz with the eigenvalue +1. 
Together with Fig. ((2}, this shows the behavior of the SNR 
between the weak-measurement regime s <^ 1 and the strong- 
measurement regime s ^ 1. 

regard the ratio 

(pz)' \/2ssina 

^Pz ^(1 + e^'^ cos a)2 + 2s cos a (cos a + e~*) 

(4.17) 
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as the SNR and we consider the optimization of this SNR. 

In Fig. m the behavior of the SNR (|4.17p is shown as a 
function of the couphng parameter s and the pre-selection 
angle a in two different ranges of s. We can see that 
there is the ridge of the SNR from the weak-measurement 
regime {s,a) — (0,7r) to the strong-measurement regime 
(s,a) = (-t-cx), 7r/2). In the strong-measurement regime, 
the fluctuation Apz in af ter th e post-selection behaves 
as ^ g/\/2s — l/{2y/{z'^)), i.e., Ap^ approach to 
zero in the limit s — > oo, while the signal (pz)' ~ 5 in 
this strong-measurement regime. Then, the SNR has the 
maximum at a = 7r/2 in the strong- measurement regime 
s ^ 1. As in Fig. [3l Fig. |4] shows the behavior of the 
SNR between the weak-measurement regime s ^ 1 and 
the strong-measurement regime s 3> 1. 

In the both of the weak measurement regime s ^ 1 
and the strong measurement regime s » 1, Fig. 0] implies 
that, for a fixed coupling parameter s, there is an optimal 
pre-selection angle a which maximize the SNR. To seek 
this optimal pre-selection angle, we consider the equation 
d{{pz)' / Apz)/da — 0. This equation yields 



cos^ a 



„coshs -I- se^ 

2 cos a 

l + s 



1 = 0. 



(4.18) 



Taking care of cos a < 1 , we easily see that the solution 
to the optimal SNR equation (|4.18|) is given by 



cos a 



cosh s + se^ 
7+1 
cos aopt(s)- 




cosh s + se^ 
7+1 



(4.19) 



This solution Q!opt('*) the pre-selection angle a which 
optimizes the SNR (|4.17|) and represents the line on the 
(s, Q!)-plane. We call this line as the optinial-SNR line. 
On this line, we can evaluate the optimally pre-selected 
SNR as 



Apz 



(s,a) 



opt 



Ap, 



■(s,aopt). (4.20) 



The optimal SNR line (|4.19p on (s, a)-plane and the op- 
timally pre-selected SNR is shown in Fig. [5] 

As shown in Fig. [5l in the strong-measurement regime 
s ^ 1, this SNR increases due to the fact that pre- 
selected state a = tt/2 is the eigenstate | fz) of the op- 
erator . On the other hand, in the weak-measurement 
regime s ^ 1, the SNR cannot be larger than that in the 
strong-measurement regime but has the minimum value 
on the optimal SNR line. Actually, for s ^ 1, the asymp- 
totic expansion of Eq. (|4.20p yields 



Apz 




0{s^), (4.21) 



which is larger than J 2/V3 1.0746 




FIG. 5: The optimal SNR line (upper panel) (|4?T9|) on the 
(s, ci)-plane and the optimally pre-selected SNR (|4.2U|I (lower 
panel) are shown. The optimally pre-selected SNR is a mono- 
tonically increasing function of s. In the strong-measurement 
regime s ^ 1, this SNR increases due to the fact that pre- 
selected state a — Tr/2 is the eigenstate | tz) of the operator 
CTz. In the weak-measurement regime s <^ 1, the SNR cannot 
be larger than that in the strong-measurement regime but has 
the minimum value (p^)' /Ap^ ~ 1.07046. 



V. APPLICATION TO THE SIMPLIFIED DSJH 
SETUP 

In this section, we apply our formulae, which are sum- 
marized in Sec. IIII A[ to the simplified DSJH [6] setup. 
We discuss the optimization of the expectation value of 
transverse defiections of an optical beam in Sec. IV Gl and 
the optimization of the SNR in Sec. IV Dl 



A. Simplified setup of experiment 

The simplified version of the DSJH experiment is the 
measurement of the tiny tilt of the piezo-driven mirror 
in a Sagnac interferometer @. In Ref. [1], they use the 
which-path information of a photon in the Sagnac inter- 
ferometer, which is represented by the photon states | O) 
and I O). Here, | O) (| C)) is the state of a photon which 
propagates along the clockwise (counter-clockwise) direc- 
tion in the Sagnac interferometer. As the pre-selected 
state of a photon, they choose 



^ (ie^'"/^] O) + e-''^/^\ 0)) , (5.1) 



V2 
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where </) is the phase difference of the states | 0) and | O) 
introduced by a Soleil-Babinet compensator. 

The weak interaction in the weak measurement is de- 
scribed by the interaction Hamihonian 



v. = kx® A5{t - to), 



(5.2) 



where k is the momentum shift of the hght path by the 
tilt of the piezo-driven mirror and x represents the shift 
of the hght image at the dark port of the interferometer. 
The quantum operator A in Eq. (15. 2p is given by 

(5.3) 



A = |0)(0|-|0)(0|. 



We note that the operator A satisfy the property A^ = 1. 
Then, we may apply our formulae given in Sec. IIIIl 

As the post-selection of a photon state, we choose the 
dark-port in the Sagnac interferometer 



IV'/)-^(|0) + ^|0)), 

and the weak value in this setup is given by 
(V-zlAIV. 



(A),„ = 



= —I cot — . 
2 



(5.4) 



(5.5) 



where (j) is the pre-selection angle in Eq. (|5.ip . We note 
that this weak value (|5.5p is pure imaginary. 

Comparing Eq. (|5.2p with Eq. (j2.ip , we find the corre- 
spondence of variables as 



X, -p, k, 



(5.6) 



where new variables x and p satisfy the commutation 
relation [x,p] = i. 

Although Dixon et al. modified the beam radius of 
the laser by lenses in Ref. we do not take account of 
the effect of this modification. This modification is not 
essential to the basic mechanism of the weak measure- 
ment. This is the reason why we call the "simplified" 
DSJH setup in this paper. 



B. All-order expectation values and probability 
distribution 



Here, we apply the formulae summarized in Sec. IIII Al 
to the above simplified DSJH setup. 

The normalization factor Z [Eq. p.3p ] , is given by 



^/{{x- (x)')^)' and Ap 



Fluctuations Aa; : 
V {{p ^ {p)')'^)' in X and p are given by 

s^e"" (cos (f) — e"") 



fc^(Ax 
1 



{Apf 



s 

2 (1 — e^'* cos ( 
1 1 

2s 1 — e^* cos (h 



\2 ' 



(5.10) 
(5.11) 



As in the case of AAV setup, the first term in Eq. (I5.10p 
[Eq. (|5.11|) ] shows the initial variance in x (in p). The 
remaining terms in Eqs. (jS.lOl) and (j5.1ip are due to the 
pre-selection, weak interaction, and the post-selection. 

In Fig. ini the expectation value —k{x)' of Eq. 
is shown as a function of the coupling parameter s 
2fc^(x^)' and the pre-selection angle (j). 



0.4 



0.4 < 

A 




FIG. 6: —k{x)' [Eq. (|5.9p ] is shown as a function of the 
coupling s — 2k^{x'^) and the pre-selection angle (j) in two 
ranges of s. (a) s G (0,10), (b) s G (0,2). For large s, the 
expectation value —k{x)' decays exponentially. From these 
figures, we can see that, for a given coupling parameter s, 
there is a pre-selection angle (p such that the expectation value 
of X is maximized. 



In the simplified DSJH setup, the probability density 
in a;-space is obtained from Eq. (13. 9p as 



1 — cos(2fcx — 4>) 
1 — e~'^ cos 6 



{x\pd\x), (5.12) 



where {x\pd\x) is the initial probability density in x- 
space: 



exp 



(5.13) 



1 — cos ( 



(5.7) 



where s is the coupling parameter [see Eq. (13. 4p ] defined 

by 



■.^2k^(x^) 



(5.8) 



The expectation values of x and p after the post-selection 
are given by 



k{xy 



se sm (p 
. — cos ( 



-, {pY^o. 



(5.9) 



C. Expectation value optimization 

Here, we consider the optimization of the expectation 
value amplification in the simplified DSJH setup. From 
Fig. m we can see that the expectation value (j5.9p expo- 
nentially decays in the strong measurement regime s 3> 1 
(Fig. EKa)). Furthermore, Fig. (B^b) also shows that, for 
a given coupling parameter s, there is a pre-selection an- 
gle 4> such that the expectation value of x is maximized. 
This is the optimal expectation value of —k{x)' for a fixed 
coupling parameter s. To seek this optimal expectation 
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value, we consider the equation d{—k{x)') /dcj) ■ 
yields the equation 



0, which expectation- value line (j5.14p is given by 



(5.14) 



This is the equation for the optimal expectation- 
value line on (0, s)-plane. On this optimal line, the 
expectation-value of x and the fluctuation Ax in x are 
given by 



k{xy 

kAx 



se 



Vl - e- 
s 



: fc(x)^pt, (5.15) 



(5.16) 



We note that the variance [Eq. (|5.16p ] in x after the opti- 
mal post-selection coincides with that of the initial state 
of the detector. The optimal expectation- value line (j5.14p 
and the expectation value of —k{x)' (|5.15p on this opti- 
mal line is shown in Fig. [T] 




FIG. 7: The optimal expectation-value line (I5.14p (upper 
panel) and the expectation value of —k{x)' (|5.15p (lower 
panel) on this optimal line is shown as a function of the cou- 
pling parameter s. In the limit s — >■ oo, the optimal pre- 
selection angle for expectation value approaches to </> — > 7r/2. 
Further, on the optimal expectation-value line, the expecta- 
tion value ~k{x)' has the maximum value at s ~ 0.8. 



{x\p'd\x) 



1 — e * cos(2/ca;) — \/l — e sin(2fca;) 



1 -e-^" 
exp 



(5.17) 



which is shown in Fig. [Sj The s = 0.1 case, which 
corresponds to the weak measurement on the optimal 
expectation- value line (I5.14p . shows that the peak of the 
probability distribution slightly deviates from the linear 
result sSA^ given in Ref. @, and that the probability 
density after the post-selection is slightly different from 
the Gaussian distribution. When the coupling parameter 
s is large, many peaks appear in the probability density 
in x-space and the expectation value (x)' after the post- 
selection approaches to zero due to the contribution of 
these many peaks. 




FIG. 8: [Online Color] Probability distribution functions 
(|5.17p in a;-space on the optimal expectation- value line (|5.14p 
with some coupling parameters s are shown by the thick line 
(red line). We also plot the initial probability distribution 
{x\pd)/k in Eq. (|5.13[l with the same coupling parameter s 
by the thin lines (blue lines). The s = 0.1 case corresponds 
to the weak measurement on the optimal expectation-value 
line (|5.14p . This shows that the peak of the probability dis- 
tribution slightly deviates from the linear result sSsA^,. The 
s — 1.0 case is still essentially same as s = 0.1 case. The max- 
imal expectation value (s = Sm ~ 0.8) is obtained around this 
parameter. The s = 10 and s = 1000 cases correspond to the 
strong measurement case. 

In the limit s — >■ oo, the optimal expectation- value 
line approaches to — )■ tt/2. Further, we have to note 
that, on the optimal expect at ion- value line, the expec- 
tation value —k{x)' has the maximum value at s 0.8. 
This is maximal value of —k{x)' on whole (0, s)-plane. 
To seek this maximum point, we consider the equation 
d{k{x)'^^^) / ds — 0. The solution s = Sm to this equation 

is derived from the equation 



From Eq. (|5.12p , the probability density on the optimal 



1 



-2s„ 



0. 



(5.18) 



12 



The numerical value of Sm is Sm — 0.79681. 
the expectation value satisfy the inequality 



Therefore, 



Hx)' < -k{xy. 



opt\s=s 



0.402371. 



(5.19) 



At this maximum point, the optimal pre-selection an- 
gle (/>m is determined by cos^m = e^*", which yields 
(j)m — 1.103 rad ~ 63.2°. We also note that the expec- 
tation value (x)' at the maximum point s — Sm itself is 
proportional to k^^. Therefore, we can obtain the large 
expectation value (x)' if we have a small coupling con- 



stant k in the interaction Hamiltonian (|5.2|) . 

If we evaluate the amplification factor A by follow- 
ing to the discussion by Dixon et al. Q, the amplifica- 
tion factor is given by ^ = \{x)'\/S. Here, 6 = klmd/ko 
is the unamplified deflection without the interferometer. 
The unamplified defiection in their experiment is (5 ~ 3 
fim. On the other hand, from Eq. (|5.19p . we obtain 
l(^)optls=s„, I — 0.4/fc at the maximum point s = s,„. 
Since fc ~ 2 x 10~^ ^ui"^ in their experiment, the maxi- 
mal amplification is estimated as Amax ^ 600. 

However, since the definition of s is given by Eq. (15. 8p 
and the fluctuation (x^) of the initial variance is re- 



garded as the beam radius, Sr, 



0.8 corresponds to 



x^) ~ 0.63 X (l/k). The optimal beam radius in their 

3 cm from Eq. (|5?T9l) . On 



Vi 

setup is given by -\/(a 
the other hand, the maximum expectation value ~{x)' ~ 
0.4 X (1/fc) ~ 2 cm. For one-photon case, the SNR at 
the maximal point s — Sm is given by |(a;)'|/Aa; ^ 0.6, 
which is independent of the coupling constant k in the 
interaction Hamiltonian. For this reason, in Sec. IV Dl we 
consider the optimization of this SNR. 



D. SNR optimization 

As in the case of the AAV setup, we consider the op- 
timization of the SNR. From Eqs. ([5J|) and (f5T0|) . the 
SNR is given by 



2se * sin < 



(1 — e ^coscf)) +2se *(cos0— e 



(5.20) 

The behavior of this SNR on (0, s)-plane is shown in 
Fig. m which indicates that the SNR (|5.20|) in the sim- 
plified DSJH setup is maximized only in the weak mea- 
surement regime s < 1. 

To carry out the optimization of the SNR (|5.20l) . we 
consider the equation d{\{x)'\/ Ax) /d(f> = 0, which yields 



cos 



' + 2 



(cosh s 



■ C0S( 



1 = 0. 



(5.21) 




v., A 



FIG. 9: The signal to noise ratio (SNR) ((5^ in simplified 
DSJH setup are shown as a function of the coupling parameter 
s and the pre-selection angle 0. This shows that the SNR in 
this setup have the peak only in the weak measurement regime 
s < 1. 



The solution (j) = 0opt Eq. (j5.2ip is given by 



cos<?!)opt(s) := 



cosh s — se 
s - 1 



y^'(cosl7; 



')2-(s-l)2 



(5.22) 



Equation (|5.22l) describes the optimal line on {(j), s)-plane. 
On this optimal line, the optimal SNR is given by 



l(=^>'l 



Ax 



\{^y\ 



opt 



Ax 



(5.23) 



Popt 



Along the optimal line [Eq. (|5.22p ]. the optimized SNR 
is shown in Fig. [TOl The optimal SNR (j5.23p is a mono- 
tonically decreasing function of s. Further, only in the 
region s < 0.15, the optimal SNR (|5.23p can be larger 
than unity. Actually, for s ^ 1, the asymptotic behavior 
of the optimal SNR (|5.23l) is given by 



Ax 



opt 



< 




(5.24) 



where ^2/V3 1.0746. Thus, we have shown that the 
upper limit of the SNR in the simplified DSJH setup for 
the single photon case is of the order of unity. 
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FIG. 10: The optimal-SNR pre-selection angle </>opt [the so- 
lution to Eq. (|5.22|) ] (upper panel) and the optimized SNR 
(|5.23p (lower panel) in simplified DSJH setup are shown as 
a function of the coupling parameter s. This shows that the 
SNR in this setup have the peak only in the weak measure- 
ment regime s < 1. The maximum SNR is of the order of 
unity in the single photon case. 



VI. SUMMARY 

In summary, after reviewing the formulation by Wu 
and Li [l3l | , we derived some formulae for the weak mea- 
surement of the operator A which satisfies the property 
A? = 1 through their formulation. We have to emphasize 
that our formulae are based on the exact evaluation of 
the formulation of Wu and Li. In the derivation of these 
formulae, we assume that the initial state of the detector 
is zero mean- value Gaussian. We note that we do not use 
any additional condition to derive these formulae. Our 
formulae are valid not only in the weak measurement 
regime but also in the strong measurement regime. Due 
to this fact, we could clarify the connection between the 
strong measurement regime and the weak measurement 
regime. 

We applied our formulae to two experimental setups. 
One is the experiment of the weak measurement using 
spin-1/2 particles, which was proposed in AAV original 
paper of the weak measurement Jj . The other is the sim- 
plified version of the optical experiment in the Sagnac in- 
terferometer (simplified DSJH setup) by Dixon et al. Q. 
These two experimental setups are typical experiments 
of the weak measurements. The weak value is real in the 
AAV setup, while it is pure imaginary in the simplified 
DSJH setup. In these setups, we have two control pa- 



rameters. One is the pre-selection angles in these exper- 
iment and the other is the coupling parameter s defined 
by Eq. (|3.4p . We discussed the behavior of the expecta- 
tion values of the detector variables in the whole range 
of these two parameters. 

In both setups of AAV and DSJH, we found that for a 
fixed coupling parameter s, there exits the pre-selection 
(or post-selection) which maximize the expectation val- 
ues of variables for the detector or the SNR. The pre- 
cise estimation of this optimal pre-selection (or post- 
selection) is possible through the exact expression of the 
expectation values summarized in this paper. This is the 
main results of this paper. The existence of this optimal 
pre-selection (or post-selection) comes from the fact that 
we specify the subensemble of system through the post- 
selection in weak measurements. Since the post-selection 
is the restriction of the system ensemble, the density ma- 
trix of the detector after the postselection is renormalized 
by this restriction. This is the essential reason of the ap- 
pearance of the normalization factor Z in Eq. (|2.9p . The 
behavior of the normalization factor Z leads the existence 
of this optimal pre-selection (or post-selection). 

Furthermore, we showed that the optimized SNR is or- 
der of unity in the weak measurement regime for the sin- 
gle particle (or photon) experiment in both experimental 
setups. To improve this SNR, we have to consider the 
large N ensemble of particles (or photons). Due to this 
large N ensemble, the SNR is improved b y th e factor V7V 
as proposed by Aharonov and Vaidman [lO|. In particu- 
lar, the photon number is very large in the experiments 
using the laser beam (for example, the experiment by 
Dixon et al. @). For this reason, the large SNR should 
be obtained in the actual experiments. 

Finally, we have to emphasize that many other exper- 
iments are also categorized into the case of = 1 and 
the initial Gaussian state of the detector. For example, 
the experiment by linuma et al. [8] corresponds to the 
experiment to measure the operator A, which satisfies 
A^ = 1, with a real weak value. They experimentally 
confirmed the formula p.Sp . The experiment by Hosten 
and Kwiat 0] corresponds to the experiment of the oper- 
ator A, which satisfies A^ = 1, with a weak value of pure 
imaginary. Thus, we may say that there are many exper- 
iments to which our formulae are applicable. Of course, 
in some actual experiments, there are some complexity 
which we did not take into account in this paper. For ex- 
ample, the modification of the beam radius by lenses in 
DSJH experiment [6] was not included in our treatment. 
Furthermore, we might have to care about the validity 
of the von-Neumann interaction model (|4.2|) in real ex- 
periments. Although there are some issues to be taken 
into account when we apply our arguments to specific ex- 
periments, we expect that our exact expressions of some 
expectation values in a weak measurement will be useful 
to understand experimental results or to propose some 
new experimental setups. 
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Appendix A: Derivations of formulae 

In this appendix, we show the derivations of the for- 
mulae summarized in Sec. IIIII Our derivation use the 
Wu-Li formulation reviewed in Sec. |lll Since their 
formulation requires the separate treatments according 



to the norm of the pre- and post-selected states, we first 
consider, in Sec. lA 1[ the non-orthogonal case in which 
the pre- and post-selection is not orthogonal. Then, in 
Sec. lA 21 we consider the case where the pre- and post- 
selected states is orthogonal. 



1. Non-orthogonal case 

When the norm KV'/IV'j)! is non- vanishing, the density 
matrix after the post-selection is given by Eq. (|2.10p . 
Only through the property = 1, the normalization 
factor Z and the density matrix p^, which are given by 
Eqs. (12. 9p and (|2.10l) . are reduced to the following series 



^ (2n) 



2n-|-l 



n=0 

I / Jln-lk „ 'Ik 



- (2,1 + 1)! 



-(p2"+i)3(A)^, 



(Al) 



ri-l 



^ (2n) 

n=l ' 



\k=0 



k=0 



+ 00 



eLtId! ((a)«.e2„+iq 



n 



n=0 



■2kP 



„2n+l-2fc „ 2k 



PdP - {A.}^, ^2n+lt'2k+lP PdP 



k=0 



k=0 

n2 o n A rt^ 



(A2) 



From this density matrix \K2\ . we can evaluate the expectation values of p, g, p , and q after the post-selection as 
follows 



^ (2n) 

71—1 ^ ^ 



(2n-f 1) 



(A3) 



Z{qy = (g)+5(5R(A)„-l-5(A)„((gp + M))) + x*S(A)^E^77e^^('^^''"^'+P'"^'9) 



(2?i + l) 



--(i-KA) 



n=l ^ ' 



= (/) + i(l-|(A)J^) E 

\n— 1 



(2n)! 



2n 



' -t-oo 



E 



(-*25) 



2n-|-l 



(2n-|-l)! 

n— ^ ^ 



2n+3\ 



(A4) 
(A5) 



= (q2) + ,^ (25R(A)^(<z) +3(A),„(<z2p + p52^) ^ |. (2|(A).„|2 - (l - |(A),„|2) (gV +^'9')) 



4 i-i(A).nE 



+ 00 



n^2 



^^1^ ((gV" +P^>^> + ^(2n)(2n - l)(p^"-^) 



+ 00 



;<,V"+' + i(2r!+ l)(2n)(p- 



2n-l\ 



2---'-A^ (271 + 1)! 

n— 1 ^ ' 

where l^^y \— Tr^/ Further, the probability densities in j>-space and in g-space are given by 



(A6) 



1 



zWd\p) = 

mp'M) = I(g|0>|'+I + J7 



1 + - ( 1 - |(A)^n (cos(25p) - 1) + 3(A)^ sin(2gp) 



{p\Pd\p), 



(A7) 
(A8) 
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where X and J are defined by 

,,2n f " / a2n~2k \ / a2k 



+ g2n J " / g2n~2k \ / g2k \ ' 



/ C)2ji-2fe-l \ / a2fc+l 



a2ji-2fe-l \ / a2 _ , _ 



fe=0 



Ml' 



a2n-2fe \ / g2k+l 

(A>;E2„+lC2.+l (^^('Zl'^)) (MO) 



fc=0 



When the initial state of the detector is zero-mean 
value Gaussian p.ip . the moments of p are given by 
Eqs. p. 21) . In this case, the normalization Z [Eq. (jAip ] 
is given by Eq. p.3p . Similar calculations with the prop- 
erties 



for n > 1, (All) 
2n- 1 (2n- 1)!! 



n-l 



2 (2a) 
0, for n > 1 



(A12) 
(A13) 



of the zero mean-value Gaussian state lead to the ex- 
pectation values of p and q [Eq. (|3.5p and p.6p ] after 
the post-selection and the variances in p and q [Eq. (|3.7p 
and p.Sp ] after the post-selection. We also note that 
the derivation of the probability density p.9p in p-space 
is straight forward, while the derivation of Eq. p. 101) is 
non-trivial. Therefore, we only explain the derivation 
Eq. (IHTTUl) . 



The initial state {q\(j>) of the detector is derived from 
the Fourier transformation of Eq. p.ip : 



(91'^) = — 



1/4 



exp 



1(1 

2 u 



(A14) 



From the definition of the Hermite polynomial [1| 



i/„(x) := (-l)"e-V2 A_(e-V2)^ 
we easily obtain 



(A15) 



^{q\<P) = (-5)-"5"/2i/„ (^^) (A16) 



This formula (|A16p is used to evaluate the derivative of 
the initial wave function (|A14p in Eqs. ([M]) and (|A10p . 

To evaluate through Eq. (|A8p . we first consider 

the second term I in Eq. (|A8P : 



-t-oo 



I = 



2fc 



(^) i?2. ^ 



n-l 



|(A)u,|^ E 2nC'2fe+l-ff2n-2fc-l H2k+1 ^ ^ 



k=0 



V 5 



(A17) 



Here, we note that the Hermite polynomial (jAlsp is Further, we also note the sum rule of the Hermite poly- 
an even function of x if the index n is even and an odd nomial |16| : 



function of x if the index n is odd llal. i.e. 



H2,,{~X) ^ H2n{x), (A18) 
H2n+l{-x) ^ -H2n+lix), (A19) 

H2k{0) = (-l)'^(2fc - 1)!!, H2k+ii0) = 0.(A20) 



1 " 



r=0 



(A21) 
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Through these formulae (IA18p - (jA2ip . we easily obtain 



X] 2kC2rH2k-2r{V2x)H2riV2x) 

= 2^-i(i/2fe(2a;) + (-l)'=(2fc-l)!!), (A22) 

2kC2r+lH2k-2r-l{V2x)H2r+liV2x) 

= 2'^-' {H2k{2x) - {-l)H2k - ly.l) , (A23) 

k 

2fc+lC'2ri?2fc+l-2r(A/2a;)ff2r(%/2a;) 

= 2^-i/2^2fe+i(2x), (A24) 

y^2fc+lC'2r+lg2fc-2r(V2x)iJ2r+l(V2a;) 

= 2^-i/2i/2fe+i(2x). (A25) 



Through the formulae (|X22]) ~ (|M5l) . Eq. (|H7)) is given 

by 



^1(910)1' 



n=0 ^ ' 



l-\{\)S\e-' 



' +CXD 



2sg 



(A26) 



Here, we note the formulae [la|: 



2 °° +2n+l 

sinh(te) = e* ^'E^2»+i(^) (2„+iV ' (^^7) 



cosh(te) 



^if2n(a 



(2r 

2n 



n=0 



(2n) 



(A28) 



Through the formula (|A28I) . 



^ = 1(910)1^ 



cosh 



2 / £9 
5 



(A),„rsinh" 



a 



(A29) 



Using Eq. (IA27[) . the similar evaluation of the final term 
J in Eq. yields 



J = |(9|0)|'3?(A).„e-^sinh ( ^ ) . (A30) 



Through Eqs. (IA29|) and (|A30|) . we can evaluate 
Eq. (|A8I) and the probability density in g-space is given 
by Eq. (I3.10p with the initial probability density in q- 
space p.l2p . 

2. Orthogonal case 

In the orthogonal weak measurements (V'/IV'i) = Oi the 
orthogonal weak value is trivial as shown in Eq. p. 131) . 
Through the orthogonal weak value (I3.13p . the normal- 
ization constant Zq defined by Eq. (I2.14p is given by 

±^ (-402)" /„2n+2x 

Zo = l + 2y ; ^ \, ' . (A31) 

^^(2n + 2)! ^ ' 

The density matrix of the detector after the post- 
selection is given by Eq. (13.151) . 

The expectation value of p, q, , and q^ after the 
post-selection are evaluated as 



z,{p'){py = z,{p')T,{pp',) 

n—1 ^ ' 

Z,{p'){qy = Z,{p')Tr{qp',) 

+ A.(2n + 2)! 

Z,{p'}{py = Z„(p2)Tr(p2p',) 



(A32) 



(A33) 



+ 00 , . 2\n 

- (P'>+2E7E^(P'"^')' (A34) 



^,(2n + 2)! 
Z,{p'){qy = Z„(p2)Tr[gVy 

-1 I 1 / 2 2 I 2 2\ 
= 1 + 2^ P +P Q ) 



_ 1 ^ ( V)" ('2(o2b2»+2 + »2n-H2 2n 

n—1 ^ ^ 

+ (2n + 2)(2n+l)(p2")) . (A35) 



From the density matrix p.l5p . we can directly obtain 
the probability density in p-space as 

{P\P» = ^^J^ (1 - cos(2<?p)) {p\pd\p)- (A36) 
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On the other hand, we also obtain the probabihty density in g-space as 



+00 



Zo{p^){q\p'd\q) ^ {q\ppd 



'dP\q) + (2n + 2)! ^ 2n+2C2/c+ 



(A37) 



since we choose the initial state of the detector as a pure state pd — \(I)){4'\j Eq. (|A37p yields 



dq 



n=l ^ ' fe=0 



^2(ji-fc) + l 
5^2(n-fc) + l 



ak.+i 



2k+l 



(A38) 



When the initial state of the detector is Gaussian p. II) . 
we use Eqs. ((3?2]) and (|A11I) - (|A13|) . Then, the expecta- 
tion values of p and q after the post-selection are trivial 
as shown in Eqs. p.l6p and p.l7p . Since the expecta- 
tion values of p and q after the post-selection vanish, the 
expectation values of p^ and q^ themselves represent the 
variances in p and q after the post-selection. Then we ob- 
tain Eqs. (|3.18p and (I3.19p . Furthermore, the probability 



density (|X36)) in p-space trivially yields Eq. (|3.20p . How- 
ever, the expression of the probability density p.2ip in q- 
space requires the non-trivial derivation from Eq. (jA38p . 
Therefore, we briefly explain the derivation of Eq. (I3.2ip 
below. 

Since our initial state of the detector is a zero mean- 
value Gaussian (jA14p . we also apply the formula (|A16I) . 
Substituting Eq. (|JT6)) into Eq. (|X38l) . we obtain 



Z,{p'){q\p'd\q) = \{qm''-^ 



Through formulae (jA28| and (|A23[) , the calculations sim- ilar to the derivation p.lOp yields p.2ip . 
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